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1.  Introduction 

This  report  is  the  final  report  under  this  project  that  was  started  on  September  11,  1992. 
The  research  has  been  carried  out  under  Grant  F19628-92-K-(X)26  from  the  Phillips  Laboratory 
at  the  Hanscom  AFB,  MA.  The  funding  was  awarded  in  response  to  a  proposal  submitted  to  the 
Lab  on  May  21,  1992.  The  original  ending  date  of  the  grant  was  8  August  1993.  It  was  later 
extended  to  1 1  December  1993. 

There  were  two  main  areas  of  study  that  were  outlined  in  the  proposal.  These  were: 

A.  Height  Bias  Model  for  Conversion  of  GPS/Geoid  Information  to  Local  Vertical  Datum 

1 .  Select  two  geographic  regions  in  the  U.S.  and/or  Canada  for  study.  Acquire 
needed  data  from  appropriate  organizations. 

2.  Develop  height  bias  model  after  examining  height  inconsistencies  assuming  all 
data  is  correct 

3.  Carry  out  estimation  of  height  bias  model  for  two  cases:  one  using  absolute 
height  determinations  and  another  using  relative  height  determinations. 

4.  Discuss  applications  of  modeling  bias  function  with  dependency  on  data 
distribution. 

5.  Define  recommended  procedures  for  bias  modeling  under  different  data  scenarios. 

B.  World  Height  System  Definition  Through  Space  Positioning 

1.  Extend  global  simulation  studies  to  incorporate  local  areas  in  which  no  global  network 
stations  are  available. 

2.  Implement  a  test  calculation  to  provide  a  vertical  datum  connection  between  two  well 
defined  vertical  datums.  We  would  propose  that  the  two  systems  be  UELN73  and 
NAVD88. 

3.  Develop  recommendations  for  the  acquisitions  and  use  of  data  for  the  determination  of 
a  World  Height  System. 

During  this  project  time  period,  substantial  research  has  been  carried  out  for  the  two  main 
areas  of  study  noted  above.  In  addition,  a  study  was  carried  out  that  used  a  large  station  position 
data  set  to  attempt  a  vertical  datum  analysis.  The  results  of  these  studies  will  be  leponed  in  this 
final  technical  report  of  the  project. 

2.  Research  Activities  and  Findings 

The  research  under  this  project  was  carried  out  by  Richard  H.  Rapp  and  Nagarajan 
Balasubramania.  Some  of  this  research  followed  ideas  developed  in  the  report  by  Rapp  and 
Balasubramania  (1992).  For  convenience,  this  final  repon  has  been  divided  into  sections 
prepared  by  each  author.  i 

2.1  Research  Activities  of  N.  Balasubramania 

In  this  project  research  has  been  carried  out  in  the  two  areas  described  in  Section  1.  The 
research  in  the  topics  of  Section  A  were  previously  described  in  Status  Report  No.  1  of  this 
nrniect  which  was  prepared  in  Januarv  1993.  For  completeness,  the  material  presented  in  the 
previous  report  will  also  be  presented  here.  The  research  under  Item  B  has  not  been  previously 
reponed  so  that  the  information  to  be  given  here  is  entirely  new. 
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2.1.1  Studies  in  Height  Bias  Model  Conversion 

The  initial  studies  in  the  project  were  related  to  the  development  of  a  height  bias  model 
using  GPS  heights,  leveling  information,  and  precise  gcoid  undulation  information.  The  research 
was  concentrated  mainly  on  the  development  of  an  height  bias  model  that  can  facilitate  the 
conversion  of  geoid/GPS  information  to  local  venical  datum.  Since  the  regional  vertical  datums 
that  are  realized  and  used  today  do  not  liave  the  ideal  geoid  as  reference  surface,  the  orthometric 
heights  computed  on  a  regional  vertical  datum  will  be  different  from  the  orthometric  height  that 
is  computed  with  respect  to  an  ideal  geoid.  The  difference  may  simply  be  a  bias  or  it  may  be  a 
difference  driven  by  a  number  of  factors  including  distortion  in  a  local  venical  network  (Rapp, 
1992).  Assuming  a  linear  relation  between  the  onhometric  height  computed  in  a  local  venical 
datum  Hd  and  the  orthometric  height  Hi  computed  using  the  relation, 

Hi  =  h-N  (1) 

where  h  is  the  ellipsoidal  height  of  the  station  measured  along  the  ellipsoidal  normal  passing 
through  the  station  and  N  the  geoidal  undulation  which  is  the  separation  between  the  defined 
ellipsoid  and  the  ideal  geoid,  we  can  write 

Hd  =  Hi  +  c(<t),  X.)  (2) 

where  c(<|>,  X)  can  be  considered  as  a  correction  term,  dependent  on  position  in  the  datum,  to 
correct  orthometric  heights  referred  to  an  ideal  reference  surface  to  the  regional  datum 
orthometric  height.  If  sufficient  points  are  available,  the  value  of  ‘c’  can  be  mapped  and 
represented  in  some  functional  form  depending  on  the  size  of  the  area  and  the  behavior  of  ‘c’. 

Development  of  an  height-bias  model  or  setting  procedures  for  mapping  the  bias  function 
value  ‘c’  in  some  functional  form  requires  GPS  observations  (height  determination  specifically), 
precise  geoid  undulation  model  and  consistent  onhometric  heights  in  a  local  vertical  datum. 
Currently  such  information  is  available  in  the  United  States,  Canada  and  in  some  part  of  Europe. 
On  our  request,  the  Geodetic  Survey  of  Canada  (A.  Mainville)  provided  the  GPS  data  for  four 
traverses:  South  Alberta  (106  stations).  North  Alberta  (51  stations).  Central  Alberta  (52 
stations),  and  Great  Slave  Lake  area  (91  stations).  In  addition,  the  orthomctric  heights  for  all  the 
GPS  stations  in  the  Canadian  Geodetic  Vertical  Datum  (CGVD-28)  and  the  high  resolution  geoid 
model  GSD  91  were  obtained.  For  the  United  States,  the  National  Geodetic  Survey  (D.  Milbert) 
provided  the  data  for  GPS  traverses  in  Florida  (52  stations),  Virginia,  G105  (62  stations),  and  in 
Oregon  (44  stations).  They  also  supplied  the  orthometric  heights  for  all  the  stations  in  National 
Geodetic  Vertical  Datum  NGVD-29  and  the  high  resolution  Geoid  Height  model  GEOID90 
(Milbert,  1991b). 

Based  on  the  data  received,  efforts  to  develop  a  suitable  height  bias  model  were  made 
during  the  period  under  reporr.  The  procedure  followed  is  as  follows: 

1.  Since  the  GPS  station  coordinates  in  the  Canadian  GPS  traverses  were  given  in  the 
WGS  84  coordinate  system,  they  were  first  transformed  to  the  ITRF  90  geocentric 
coordinate  system  using  the  transformation  parameters  given  in  Boucher  and 
Altamimi  (1991),  to  make  the  origin  of  the  coordinate  system  compatible  with  the 
geopotential  model  used.  In  the  case  of  GPS  traverses  in  the  United  States,  the  station 
coordinates  were  already  in  the  ITRF  90  geocentric  system  as  reported  by  Milbert 
(1991a)  and  hence  no  transformation  was  required.  The  geoidal  undulation  at  these 
stations  were  then  interpolated  using  the  respective  geoid  models. 
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2.  The  orthometric  height  bias  ‘c’  was  then  computed  using  the  relation. 


c((t),  X,)  -  (hcps  -  Nceoid  Model)  *  Hiocal  datum  (3) 

Statistical  analysis  of  the  height  bias  values  for  US  and  Canadian  traverses  are  given 
in  Table  1.  An  example  of  an  orthometric  height  bias  table  is  given  in  Table  2  for  the 
Slave  Lake  Area. 

Table  1;  Statistics  of  the  Difference  in  Orthometric  Heights  Derived  from  GPS/GSD91  and 
CGVD28  Heights  in  case  of  Canadian  Traverses  and  from  GPS/GEOID90  and  NGVD29  Heights 
for  GPS  Traverses  in  the  United  States.  Units  are  in  cm. 


Traverse 

United  States  1 

Canada 

ijm 

j  Florida 

Ksuial 

North 
Alberta  I 

Central 

Alberta 

1  South 

1  Alberta 

||lg?ngn?cQ| 

#  of  Stations 

44 

HM 

wmm 

45  j 

■EUH 

5i 

i  91 

i  -77 

-64 

73  I 

WMSM 

. "49 . 

wmm 

[  -31 

Max.  Ht.  Bias  Value  j 

+58^ 

^  +49 

153  1 

\  +62  1 

-45 

Hcm 

+48 

Mean 

2 

3“ 

1  -5 

119  1 

1  -1  i 

-59 

i  23 

-3 

Std.  Deviation  I 

22 

39 

39 

1 

11  1 

8 

21 

'  16 

3.  Three  different  procedures  were  attempted  to  map  the  height  bias  function  value  ‘c’. 
They  were: 


•  global  surface  fitting  using  deterministic  polynomial  functions, 

•  minimum  curvature  fitting  using  cubic  spline  functions,  and 

•  height  bias  modelling  using  least  squares  collocation  technique. 

Since  the  polynomial  functions  provide  unpredictable  errors  in  the  interpolated  values  away  from 
the  stations  and  also  the  minimum  curvature  technique  do  not  provide  the  estimated  accuracy  of 
predicted  height  bias  function  values,  the  least  squares  collocation  technique  was  attempted  for 
modelling  the  height  bias  function  value  ‘c’.  Assuming  a  correlation  length  of  40  km  and  using  a 
second-order  Markov  model  covariance  function  fitted  to  local  data,  the  height  bias  value  along 
with  its  estimated  accuracy  were  predicted  for  any  other  location  in  the  area  using  least  squares 
collocation  technique.  Contour  plots  showing  predicted  height  bias  function  in  the  area  and  its 
estimated  accuracy  one  each  for  U.S.  traverse  (in  Oregon)  and  Canadian  traverse  (in  South 
Mbena)  are  shown  in  Figures  1  (a,  b)  and  2  (a,  b)  respectively  with  a  contour  interval  of  5  cm. 
From  Figures  1(b)  and  2(b),  we  see  that  the  estimated  accuracy  of  predicted  height  bias  function 
is  of  the  order  of  20  cm  if  we  have  well  distributed  GPS  stations  in  the  area,  'Hie  advantage  of 
developing  such  an  height  bias  model  is  that,  by  knowing  the  onhometric  height  with  respect  to 
an  ide^  geoid  using  Geoid/GPS  observations  using  the  relation  (1),  its  orthometric  height  in  the 
local  vertical  datum  can  bei  computed  using  the  interpolated  bias  function  ‘c’.  Also,  the  plot 
showing  the  estimated  accuracy  of  predicted  ‘c’  values,  serves  as  a  reliability  diagram  intimating 
the  user  the  accuracy  of  the  pr^cted  ‘c’  value  and  hence  the  onhometric  height  in  local  verticid 
datum. 


Considering  Figure  la,  one  sees  a  clear  large  signature  that  is  associated  with  the  higher 
elevntions  in  the  western  part  of  Oregon.  These  lareer  bias  terms  may  be  associated  with  geoid 
undulation  errors  in  Geoid  90  associated  with  the  indirect  effect. 
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Table  I  i  Computation  of  Orthometric  height  bias  function  in 
Great  Slave  Lake  Area,  Canada.  Heights  in  meters 


Lat 

(deg) 

c-(h-N)-H 

Lon  Ellip.ht  geoid.ht  (h-N)  Ortho. ht  Ht. error 

(deg)  (h)  (N)  (H)  (in  cm) 

(from  GSD91)  (in  CGVD28) 

62.475 

245.559 

183.350 

-25.905 

209.255 

209.229 

2.5 

60.038 

243.117 

279.969 

-19.840 

299.809 

299.902 

-9.2 

60.426 

243.647 

263.951 

-21.571 

285.522 

285.730 

-20.7 

60.800 

243.408 

244.827 

-21.935 

266.762 

267.028 

-26.5 

60.709 

244.995 

216.481 

-24.309 

241.290 

241.269 

2.0 

62.479 

245.273 

155.184 

-25.623 

180.807 

180.706 

10.1 

62.544 

245.010 

141.715 

-25.083 

166.798 

166.700 

9.8 

62.593 

244.813 

142.511 

-24.582 

167.093 

167.061 

3.2 

62.660 

244.705 

152.975 

-24.274 

177.249 

177.221 

2.8 

62.686 

244.535 

148.808 

-23.980 

172.788 

172.753 

3.4 

62.709 

244.411 

147.611 

-23.776 

171.387 

171.339 

4.8 

62.784 

244.007 

140.531 

-22.935 

163.466 

163.455 

1.0 

62.641 

243.743 

210.678 

-22.284 

232.962 

232.996 

-3.3 

62.556 

243.591 

233.449 

-22.033 

255.482 

255.507 

-2.5 

62.382 

243.503 

250.968 

-21.915 

272.883 

272.957 

-7.3 

62.302 

243.571 

221.520 

-22.063 

243.583 

243.699 

-11.6 

62.204 

243.683 

206.345 

-22.334 

228.679 

228.722 

-4.3 

62.109 

243.704 

191.012 

-22.422 

213.434 

213.524 

-9.0 

62.024 

243.687 

188.108 

-22.451 

210.559 

210.650 

-9.1 

61.979 

243.563 

196.285 

-22.232 

218.517 

218.608 

-9.1 

61.837 

243.329 

203.860 

-21.826 

225.686 

225.825 

-13.8 

61.781 

243.228 

201.425 

-21.666 

223.091 

223.148 

-5.7 

61.711 

243.103 

197.065 

-21.465 

218.530 

218.533 

-0.3 

61.677 

242.975 

190.166 

-21.261 

211.427 

211.571 

-14.4 

61.601 

242.859 

175.645 

-21.128 

196.773 

196.921 

-14.7 

61.514 

242.754 

154.450 

-21.025 

175.475 

175.646 

-17.1 

61.428 

242.601 

138.736 

-20.823 

159.559 

159.752 

-19.2 

61.366 

242.499 

136.108 

-20.732 

156.840 

157.021 

-18.0 

61.095 

242.500 

174.119 

-20.585 

194.704 

194.958 

-25.4 

60.870 

243.269 

206.432 

-21.749 

228.181 

228.496 

-31.4 

60.717 

243.528 

231.687 

-22.030 

253.717 

253.948 

-23.0 

60.629 

243.659 

240.712 

-22.095 

262.807 

262.961 

-15.3 

60.784 

245.342 

211.651 

-25.497 

237.148 

237.128 

2.0 

60.611 

245.510 

245.756 

-25.469 

271.225 

271.277 

-5.2 

60.536 

245.610 

222.953 

-25.550 

248.503 

248.527 

-2.4 

60.407 

245.733 

242.665 

-25.609 

268.274 

268.213 

6.1 

60.228 

246.153 

247.940 
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FIGURE  1A 

HEIGHT  BIAS  FUNCTION  WITH  THE  OREGON  DATA  SET 
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FIGURE  IB 

ACCURACY  OF  ESTIMATE  OF  PREDICTED  HEIGHT  BIAS  FUNCTION 

WITH  THE  OREGON  DATA  SET 
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FIGURE  2A 

HEIGHT  BIAS  FUNCTION  WITH  THE  SOUTH  ALBERTA  DATA  SET 
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FIGURE  2B 

ACCURACY  ESTIMATE  OF  PREDICTED  HEIGHT  BIAS  FUNCTION 

WITH  THE  SOUTH  ALBERTA  DATA  SET 
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2.1.2  Studies  in  the  World  Height  System  Definition 

Basically  four  kinds  of  data  are  essential  for  modelling  the  World  Height  System  (WHS). 
They  iu'c:  accurately  determined  free-air  gravity  anomalies,  precise  orthometric  heights  of  the 
stations  above  the  regional  vertical  reference  datums,  ellipsoidal  heights  of  the  stations  above 
adopted  goecentric  ellipsoid  and  accurate  geopotential  models.  Since  such  data  availability  over 
the  world  are  neither  uniform  nor  to  the  required  precision,  different  modelling  procedures  are 
needed  to  be  developed  to  accommodate  the  available  heterogeneous  data  to  come  up  with  an 
elegant  observation  equation  for  estimating  the  parameters  defining  a  WHS.  As  understandable, 
it  may  require  several  iterations  before  realizing  an  ideal  WHS  which  can  meet  the  present  day 
requirements  of  WHS  established  to  an  accuracy  of  ±10  cm. 

The  two  different  data  scenarios  anticipated  in  developing  the  procedure  for  setting  up  the 
observation  equations  for  defining  the  WHS  are: 

(a)  Given  the  free-anomalies  in  the  Molodensky  sense  and  normal  heights  of  the  stations,  and 

(b)  Given  the  free-air  anomalies  in  the  classical  sense  and  onhometric  heights  (or  normal 
orihometric  heights)  of  the  stations. 

Detailed  discussion  on  the  definition  and  comparison  of  free-air  anomalies  obtained  in 
Molodensky  sense  and  Classical  sense  can  be  found  in  Sideris  and  Forsberg  (1991)  and  for 
various  height  definitions  given  in  Rapp  and  Balasubramania  ( 1992). 


P 


Fig.  3  Geometric  representation  of  various  equipotential  surfaces 

In  Figure  3,  Qj  is  the  fundamental  benchmark  serving  as  the  origin  for  the  regional 
vertical  datum  'i'.  The  potential  of  the  reference  surface  for  WHS  is  Wq.  We  define 

C  p  =  Wo-W^  as  the  potential  difference  between  the  reference  equipotential  surfaces  for 

'*1 

datum  'i'  and  the  WHS  reference  surface.  Introducing  a  reference  ellipsoid  where  Uq  is  the 
potential  on  the  surface,  we  define  AWq  =  Wq  -  Uq.  The  distance  between  the  equipotential 
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surfaces  U  =  Uq  and  W  =  Wq  is  also  shown  in  the  figure,  and  depends  on  the  disturbing  potential 
at  P'  and  AWq. 

Procedure  for  setting  up  the  observation  equation  for  case  (a)  is  as  follows: 

Let  H*  be  the  normal  height  of  the  point  P  with  respect  to  the  WHS  (to  be  defined)  and 
H*  the  known  normal  height  of  the  same  point  referred  to  i^"  regional  venical  datum.  Then  from 
equation  (2-122)  in  Pavlis  (1988),  we  can  write  the  relation; 


AgM=Agii> 


W-. 


Yq, 


^dh 


Qo 


H- 


-QiO 


(4) 


where  Ag^.Agji'  are  the  gravity  anomalies  computed  with  normal  gravity  yg  determined  using 
normal  heights  H*,H*  of  the  point  P  referring  to  WHS  and  regional  vertical  datum  'i' 
respectively. 

From  Moritz  (1980,  eq.  48-1),  the  computation  of  height  anomaly  Cp  by  analytical  continuation 
to  point  level  can  be  done  using  the  formula: 


Cp  =  Co  +  4|^JJ(Ag^+gi)S(V)do 


(5) 


where 


r  SGM 
^0  =  -^ 


AWq 

Y 


Ag5(  is  the  free-air  anomaly  defined  in  Molodensky  sense  and  corrected  for  systematic 
effects  as  discussed  in  Rapp  and  Pavlis  (1990), 

gj  is  the  Molodensky  correction  term, 

R  is  the  mean  radius  of  earth,  y  the  mean  gravity,  and 
S(v)  the  Stokes'  function. 

If  Faye  anomalies  are  available  equation  (5)  can  be  replaced  with  eqn.  (24)  in  Wang  (1993)  as: 
C,  =  Co  +5|^J/(Ag&  +C)S(V»Xio-.^-I^H^  (6) 

where  Hp  is  the  normal  height  of  point  P  and 

5h2  =  0.453  -  0.018  sintft  +  0.087  cos<|)cosX.  -»•  0.204  cos(|k;osA.  (km^) 

Combining  equations  (4)  and  (6),  we  get 
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(7) 


;,=;.+j^Jj[igM"+c]s(v)do-j^iq  C<,  JS(V)do, 


•^£h  |_j*  ^Gp  g^2 

y  '  y  '  y 


where  qJ  =  — L- 
Vqo 


+fZl 

o.  V^h^o. 


^ahy 


H* 


From  equation  (2-121)  and  (2-122)  in  Heiskanen  and  Moritz  (1967)  we  can  write 


qi  =--^ 

^  R 


(8) 


From  Figure  3  we  note  that: 


h;  =  H>^^andhp  =  H*p -(-Cp  orCp  =  hp  -  H*p 


(9) 


where  hp  is  the  ellipsoidal  height  of  the  point  P  and  Cp  its  height  anomaly.  Comparing  eqns.  (7) 
and  (9)  and  rearranging  known  and  unknown  terms  we  write: 


Y  =  h-H:-Cl=-^^  + 

Y  y 


I  j-i 


47t 


JS(i|r)daj 


(10) 


where  +CjS(\|/)da--^^H*--^^^H*^  -  ^^^5h^  is  the  gravimetric  height 

anomaly  (without  term)  and  referred  to  regional  vertical  datum  'i'. 

Similarly  for  case  (b)  when  the  free-air  anomalies  in  classical  sense  are  given  along  with 
orthometric  heights  of  stations,  modelling  procedure  given  in  Rapp  and  Balasubramania  (1992) 
can  be  followed  to  come  up  with  the  observation  equation: 


Y  =  h-H.-N| 

•  Y  Y  Y^; 


i/SWdOi 


(11) 


where  N;  =;~pjJ(Ag'  -HC)S(v)do-i^^Hf 
term)  and  referred  to  regional  datum  'i'. 


is  the  gravimeuic  geoid  undulation  (without  No 


Equations  (10)  and  (1 1)  are  the  key  observation  equations  for  solving  the  unknown  parameters 
AWoUndCq^o  (for  i  =  1,  2,  ...  I).  The  parameters  can  be  obtained  using  a  least-squares 
adjustment  procedures.  The  adjustment  process  discussed  by  Rapp  and  Balasubramania  (1992, 
Section  4.4)  can  be  used  to  estimate  the  parameters  defining  the  World  Height  System. 

The  solution  of  the  parameter  vector  is  given  by 


(12) 


X  =  ( A^Z;‘ A  +  K^K)'‘  A^l;*  Y 

and  the  error  covariance  matrix  of  the  estimated  parameters  would  be 

Z^  =  (atz;‘A  +  K^K)’‘atZ;‘A(atZ;'A  +  K^K)'‘  (13) 

where  Y  is  the  nxl  observation  vector  from  equation  (10), 

A  is  the  nxm  coefficient  matrix  of  observation. 

X  is  the  mxl  parameter  vector, 

K  is  the  Ixm  coefficienrmatrix  of  the  constraints,  : 

is  the  variance-covariance  matrix  of  observations. 

Equation  (12)  and  (13)  correspond  to  equations  (4-33)  and  (4-34)  in  Rapp  and  Balasubramania 
(ibid). 

V 

Numerical  Investigation 

A  first  iteration  test  computation  to  estimate  the  parameters  defining  the  WHS  was 
carried  out  with  available  data  from  six  regional  vertical  datums.  The  different  vertical  datums 
considered  are:  Nonh  American  Vertical  Datum  88  (NAVD88),  ODN  (UK),  IGN69  (France), 

NN  (Germany),  AHD71  (Australia)  and  Scandinavian  datum.  Though  there  is  no  single  region^  • 
vertical  datum  known  as  the  Scandnavian  datum,  the  individual  datums  RH70  (Sweden),  N60  f 
(Finland)  and  NN1954,  NNN1957  (Norway)  were  combined  to  serve  as  a  single  regional  datum 
for  computational  convenience. 

The  1TRF91  coordinate  system  was  selected  as  the  consistent  coordinate  system  for 
defining  the  space  geodetic  station  coordinates  used  in  the  realization  of  WHS.  An  ideal 
ellipsoid  with  the  following  parameters  was  adopted  as  the  reference  ellipsoid  for  the  study: 

ac  =  6378136.3  m 
f  =  1/298.257222101 
GM  =  3986006  X  lO^mV^ 
a)  =  7292115  x  10-*lradS-‘ 

The  last  three  parameters  are  the  same  as  the  parameters  of  the  GRS80  reference  ellipsoid.  * 

Table  3  below  lists  the  distribution  of  space  stations  in  various  regional  vertical  datums  and 
Figure  4  shows  their  locations. 

I  'fable  3:  Distribution  of  space  stations 


■KliH 

1  No.  of  otauons  j 

1  3  1 

!  i 

flHlj 

The  accuracy  and  grid  size  of  teixestrial  gravity  data  used  varied  from  one  regional  datum  to  * 
another.  For  North  .American  datum.  3'x3'  atmospherically  connected  Helmert  anomalies  "■ 
around  each  of  the  U.S.  stations  received  from  D.  Milbert,  NGS  were  used.  As  quoted  by 
Milben  (1991)  these  3'x3'  mean  anomalies  have  a  random  error  of  ±1.5  mgal.  For  the 
Scandinavian  datum.  Faye  anomalies  received  from  Dr.  Rene  Forsberg  were  used.  For  the  other 
stations  in  England,  Germany,  France  and  Australia,  the  6'xlO'  and  2'x2'  atmospherically 
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corrected  free-air  anomalies  referring  to  OSU  GRS  (a  =  6378136.0,  f  =  1/298.257222101)  used 
by  Despotakis  (1983)  in  his  thesis  work  were  used.  These  gravity  data  sets  are  not  corrected  for 
terrain  and  as  given  in  Torgd  (1983),  they  have  a  random  error  of  ±7  mgal. 

The  elevation  of  the  space  stations  above  their  local  vertical  reference  datum  were 
obtained  from  NASA  space  geodesy  program  -  catalogue  of  site  information  (NASA  Tech  Memo 
4482).  For  the  stations  in  USA,  their  heights  referred  to  NGVD29  height  system.  Using  the 
VERTCON  software  received  from  Zilskoski  at  NGS,  the  heights  were  converted  to  NAVD88 
heights.  For  stations  in  England  and  Australia,  the  heights  given  were  normal  orthometric 
heights  and  for  stations  in  Germany,  France  and  Scandinavian  countries,  they  are  normal  heights 
referred  to  their  regional  vertical  datum.  The  global  geopotential  model  OSU91 A  to  degree  360 
was  used  as  the  reference  field  in  the  study. 

Three  different  techniques  were  used  to  compute  the  gravimetric  height  anomaly  and 
the  gravimetric  undulation  Nj  in  equations  (10)  and  (11)  respectively.  They  were: 

(i)  Modified  Stokes'  technique 

(ii)  Least  squares  collocation  technique  (LSC) 

(iii)  Interpolation  from  regional  geoid  height  models  developed  using  Fast-Fourier 
transform  technique. 

Comparison  of  gravimetric  undulations  computed  using  different  techniques  at  the  U.S.  space 
geodetic  stations  are  given  in  Table  4.  The  interpolated  geoid  height  values  from  recent 
GEOID93  Geoid  Height  model  received  from  Milbert  at  NGS,  are  also  included  for  better 
comparison.  The  differences  between  the  gravin^tric  undulation  and  geometric  undulation  are 
also  shown.  In  Table  5,  gravimetric  height  anomaly/undulation  computed  at  the  different 
European  space  stations  are  given.  Since  we  had  only  the  regional  geoid  height  model  NKG-89 
for  Scandinavian  countries,  only  for  stations  in  that  area,  the  interpolated  geoid  height  are 
included  in  Table  5.  The  gravimetric  undulation  at  the  Australian  stations  computed  using  the 
first  two  techniques  and  their  comparison  with  geometric  undulation  are  given  in  Table  6. 

From  Tables  4  and  5  we  see  that  the  rms  difference  between  the  geometric  undulation  and 
the  gravimetric  undulation  computed  using  different  techniques  are  almost  the  same  but  in  Table 
6  we  see  they  differ  by  about  60  cm.  Also  from  Table  6  we  note  the  difference  between  the 
gravimetric  and  geometric  undulations  computed  at  the  Yarragadee  station  (located  in  Western 
Australia)  and  at  Canberra  (located  in  Eastern  Australia)  differ  by  about  1 . 1  meters  though  both 
the  stations  are  considered  to  lie  in  the  same  regional  vertical  datum  (AHD71).  This  difference 
may  be  due  to  the  suggested  rise  of  mean  sea  level  up  the  East  Australian  coast  of  about  1.5  m  by 
Morgan  (1992)  author  the  same  extent  due  to  the  long  wave  length  error  in  the  global 
geopotential  model  used. 


( 
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Table  4:  Comparison  of  gravimetric  undulation  computed  using  different  techniques  at  the  U.S. 
space  geodetic  stations.  Units  are  meters. 
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Table  5:  Comparison  of  gravimetri(rhetghcanomaly  computed  using  different  techniques  at  the 
European  spacd-geoderic  stations.  *Units  are  metery. 


« 

Table  6:  Comparison  of  gravimetric  undulation  computed  usihig  different  techniques  at  the 
Australian  space  geodetic  stations.  Units  are  meters. 


Using  the  data  at  18,  s'pace  geodetic  stations  in  6  regional  vertical  datums,  and  with  7 
unknowns,  observation  equations  and  normal  equations  were  formed.  From  equations  (12)  and 
(13)  the  parameters  defining  the  WHS  and  their  accuracy  estimates  were  computed  for  different 
techniques  used  in  computing  gravimetric  height  anomaly/undulation.  The  estimated  parameters 
are  given  in  Figures  5  and  6. 
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Figure  5:  Regional  Vertical  Datum  separation  from  defined  WHS  (units  in  cm)  (from  Modified 
Stokes'  technique). 


Figure  6;  Regional  Vertical  Datum  separation  from  deHned  WHS  (units  in  cm)  (from  LSC 
technique)  ‘ 


From  Figures  S  and  6  it  is  clear  that  the  different  regional  vertical  surfaces  used  in  the 
definition  of  WHS  show  the  same  pattern  of  their  relative  positions  to  one  another  in  spite  of 
using  two  different  techniques  in  their  estimation  but  the  numerical  values  are  different.  The 
Scandinavian  datum  which  is  61  cm  above  the  defined  WHS  in  Figure  6  is  only  38  cm  above  the 
same  in  Figure  6.  This  may  be  due  to  using  limited  terrestrial  gravity  data  (y  =  1®)  for 
computing  the  gravimetric  height  anomaly  at  the  Scandinavian  stations.  Hie  Australian  height 
datum  AHD71  which  is  about  106  cm  above  the  defined  WHS  in  Figure  5  is  only  54  cm  in 
Figure  6.  The  reason  is  that  gravimetric  undulation  computed  at  the  Australian  space  station 
using  Modified  Stokes'  technique  and  LSC  technique  differs  by  about  70  cm.  This  may  be  due  to 
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error  in  the  gravity  anomaly  data  used  around  these  stations.  The  datum  differences  between 
DHNN  and  IGN69  as  intimated  by  Boucher  (private  communication)  is  about  50  cm  which 
agrees  well  with  the  differences  shown  in  Figure  5  and  6.  Except  for  the  Australian  Height 
datum,  the  differences  between  different  regional  datums  as  shown  in  Figure  5  compare  well 
with  the  estimates  given  by  Rapp  (1993). 

2.2  Research  Activities  of  R.H.  Rapp 

During  this  grant  Professor  Rapp  utilized  existing  Doppler  positioned  stations  to  study 
several  venical  datums  with  the  aid  of  gravimetric  geoid  undulations  derived  from  a  merger  of 
new  JGM-2  potential  coefficient  model  (degrees  2-70)  and  the  OSU91A  model  (degrees  71  to 
360).  This  research  led  to  a  paper  that  has  been  accepted  for  publication  in  the  Bulletin 
Geodesique  or  manuscripta  geodaetica.  This  paper  is  made  part  of  this  repon  and  starts  on  the 
following  page. 
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2.2.1  Separation  Between  Reference  Surfaces  of  Selected  Vertical  Datums 

Richard  H.  Rapp 
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November  1993 


Abstract 

This  paper  discusses  the  separation  between  the 
reference  surface  of  several  vertical  datums  and  the 
geoid.  The  data  used  includes  a  set  of  Doppler 
positioned  stations,  transformation  parameters  to 
convert  the  Doppler  positions  to  ITRF90,  and  a 
potential  coefficient  model  composed  of  the  JGM-2 
(NASA  model)  from  degree  2  to  70  plus  the  OSU9IA 
model  from  degree  71  to  360.  The  basic  method  of 
analysis  is  the  comparison  of  a  geometric  geoid 
undulation  derived  from  an  ellipsoidal  height  and  an 
orthomcu-ic  height  with  the  undulation  computed  from 
the  potential  coefficient  model  The  mean  difference 
can  imply  a  bias  of  the  datum  reference  surface  with 
respect  to  the  geoid.  Vertical  datums  in  the  following 
countries  were  considered:  England,  Germany,  United 
States,  and  Australia.  The  following  numbers  represent 
the  bias  values  of  each  datum  after  adopting  an 
equatorial  radius  of  6378136.3m:  England  (-87  cm), 
Germany  (4  cm).  United  States  (NGVD29  (-26  cm)), 
NAVD88  ^72  cm).  Australia  AHD  (mainland.  -68  cm); 
AHD  (Tasmania,  -98  cm).  A  negative  sign  indicates 
the  datum  reference  surface  is  below  the  geoid.  The  91 
cm  difference  between  the  dauiihs  in  England  and 
Germany  has  been  independently  estimated  as  80  cm. 
The  30  cm  difference  between  AHD  (mainland)  and 
AHD  (Tasmania)  has  been  independently  estimated  as 
40  cm.  These  bias  values  have  been  estimated  from 
data  where  the  geometric/  gravimetric  geoid  undulation 
difference  standard  deviation,  at  one  station,  is  typically 
±100  cm.  although  the  mean  difference  is  determined 
more  accurately. 

The  results  of  this  paper  can  be  improved  and  expanded 
with  more  accurate  geocentric  station  positions,  more 


accurate  and  consistent  heights  with  respect  to  the  local 
vertical  datum,  and  a  more  accurate  gravity  field  for  the 
Earth.  The  ideas  developed  here  provide  insight  on  the 
determination  of  a  world  height  system. 

Introduction 

This  paper  is  written  to  estimate  the  linear  separation 
between  the  origin  reference  surface  of  several  vertical 
datums.  The  process  for  doing  this  inherently  implies  a 


concept  of  a  world  height  system  as  discussed  by  Rapp 
and  Balasubramania  (1992).  This  paper  is  not  meant  to 
be  an  extensive  analysis  of  theoretical  procedures 
involved  in  vertical  datum  definition  and  connection. 
Recent  papers  in  this  area  include  those  by  Xu  and 
Rummel  (1991)  and  Heck  and  Rummel  (19^).  This 
paper  is  not  meant  to  show  results  of  the  latest 
measurement  procedures.  Examples  of  such  papers  are 
those  of  Rizos,  Coleman  and  Ananga  (1991)  and  Pan 
and  SjOberg  (1993).  This  paper  is  designed  to  examine 
some  earlier  ideas  (Rapp,  1983)  using  the  latest  gravity 
field  models  and  terrestrial  reference  frame 
transformations. 

Principles 

Let  N  be  the  geoid  undulation  referred  to  a  geocentric 
reference  frame,  an  ellipsoid  of  a  specified  flattening, 
and  an  equatorial  radius  considered  to  be  an  optimum 
estimate.  (Optimum  may  be  a  matter  of  opinion  so  this 
term  is  loosely  used  here.)  Our  calculation  of  N  will  be 
done  through  a  potential  coefficient  model  of  degree  L. 
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Given  the  fully  normalized  coefficients  Cnm.  Snm*  one 
has 

yr  n=2Uy 

i  ( cos  mX  +  S„„  sin  mX)P„„  (cos  0)  +  c 

m  sO 

where  the  usual  definitions  (Rapp  and  Pavlis,  1989) 
apply.  The  c  term  represents  the  downward 
continuation  effects  from  an  evaluation  of  the  series  at  a 
valid  point  in  space  to  a  point  on  the  geoid.  We  will 
assume  that  c  is  zero-in  this  analysis.  The  even  degree 
zonal  harmonics  in  (1)  refer  to  the  adopted  reference 
ellipsoid. 

Let  h  be  the  ellipsoidal-  height  of  a  point  computed 
from  rectangular  coordinates  in  a  geocentric  reference 
frame  in  a  true  scale  system.  Let  H  be  the  onhometric 
height  with  respect  to  the  geoid  that  is  described  by  the 
N  values  of  cq.  (1).  Then: 

h  =  H  -t-  N  (2) 

If  one  were  given  h  and  N  then  H  can  be  computed. 
The  accuracy  of  the  determination  of  this  H  will  depend 
on  the  accuracy  of  h,  N  and  the  optimum  equatorial 
radius. 

Consider  next  an  orthomeuie  hefght  (Hq)  referred  to  a 
specific  vertical  datum  D.  The  ideal  situation  is  that  the 
vertical  datum  is  defined  by  an  equipotential  surface, 
near  mean  sea  level,  at  one  specific  point.  In  practice, 
vertical  datums  may  have  bem  defined  by  fixing  mean 
sea  level  to  zero  onhometric  elevation.  In  this  case,  the 
reference  surface  and  the  H  values  have  some  distortion 
with  respect  to  the  ideal  case.  In  this  paper,  we  will 
assume  that  a  vertical  datum  is  defined  by  a  unique 
reference  surface.  This  surfacr  may  not  and  actually 
will  not  coincide  with  the  ideal  reference  surface,  the 
geoid. 

Let  B  be  the  bias  between  the  ideal  system  and  a 
specific  system  such  that: 

H  =  Hd  +  B  ‘  (3) 

where  H  is  the  orthometric  height  in  the  ideal  system 
and  Hd  is  the  corresponding  height  in  the  datum 
system.  A  positive  B  indicates  the  specific  datum 
reference  surface  is  above  the  ideal  reference  surface. 
Wc  can  estimate  B  by  substituting  (3)  into  (2): 

B  =  (h-HD)-N  '  '  (4) 

t' 

The  evaluation  of  B  should  be  done  using  numerous 
stations  connected  to  the  datum  D.  Averaging 


individual  values  can  reduce  the  noise  inherent  in  the 
esumatc  of  the  three  quantities  on  the  right  hand  side  of 
cq.  (4).  A  comparison  of  B  from  different  datums  will 
yield  information  on  the  relative  positions  of  these 
vertical  datums. 

Data 

The  stations  to  be  analyzed  for  this  paper  have  been 
positioned  through  Doppler  observation  techniques 
during  the  time  period  1971-1986.  A  description  of 
aome  of  the  processing  techniques  used  in  the 
estimation  of  the  precise  positions  of  an  analogous 
i  aBtion  set  may  be  found  in  Weigel  (1993,  p.  14).  With 
r  tBs  station  set,  "heights”  above  mean  sea  level  are 
generally  available.  These  heights  are  not  precisely 
defined  in  terms  of  datum  reference,  type  of  orthometric 
height,  or  if  in  fact  they  may  be  normal  heights. 
AHiough  definitive  results  of  this  type  of  investigation 
rc({uire  such  information,  this  paper  is  considered  a 
demonstration  study  and  we  will  assume  that  the  mean 
sea  level  heights  are  orthometric  heights  connected  to 
the  vertical  datum  of  the  region  or  country.  Ultimately, 
we  will  restrict  our  analysis  to  stations  that  can  be 
associated  with  a  specific  vertical  datum. 

,.<Bcfcae  any  analysis  is  carried  out,  the  ellipsoidal  height 
givccdn  the  original  data  record  is  transferred  to  the 
ITRI90  system  (Boucher  and  Altamimi,  1992,  Table  4 
(NWL  9D))  and  referred  to  an  ellipsoid  whose 
equatorial  radius  is  6378136.3m  and  flattening  is 
1/298257. 

Analyses  .such  as  thc.se  arc  made  more  feasible  as  the 
knowledge  of  the  Earth's  gravity  field  improves.  One 
of  the  recent  models  is  JGM-2  (Lcrch  et  al.,  1993) 
whichil  a  combination  model  complete  to  degree  70. 
This  model  has  been  augmented  by  the  OSU91 A  model 
(Rapp,  Wang,  Pavlis.  1991)  from  degree  71  to  degree 
360.  AUiough  all  tests  to  be  described  in  this  paper 
have  bev  carried  out  with  both  the  OSU9I A  and  JGM- 
2  (augmanted)  model  only  (with  one  exception)  results 
with  the  most  current  model  (JGM-2)  will  be  given. 
The  geoid  undulation  can  be  computed  from  this  data 
with  an  estimated  accuracy  of  ±37  cm  (Rapp,  Wang, 
Pavlis.  f.  64,  ibid.)  globally,  which  may  be 
considerdily  poorer  in  areas  lacking  adequate  gravity 
information. 

Results  •  Global 

The  first  computation  involved  the  examination  of  all 
Doppler  pMitions  independent  of  their  geographic 
location.  (JUulation  comparisons  were  made  with  both 
the-  OSU9A  and  JGM-2  (augmented)  potential 
coefficient  models.  In  this  analysis,  stations  with 
residuals  exceeding  2m  in  absolute  value  were  rejected 
from,  the  statistical  computations.  This  2m  criteria  is 
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tighter  than  the  4in  criteria  used  in  past  studies,  but  it  is 
consistent  with  the  criteria  adopted  for  the  regional 
studies  to  be  reported  shortly.  The  2in  criteria  is 
basically  twice  the  standard  deviation  of  the 
geometric/gravimetric  difference  and  thus  is  a  2  sigma 
criteria.  The  statistics  on  the  differences  are  given  in 
Table  1  where  the  mean  difference  is  defined  through 
the  geometric  undulation  minus  the  potential  coefficient 
implied  undulation.  The  total  number  of  stations 
considered  was  2033.  The  results  indicate  little 
preference  for  one  geopotential  model  although  the  use 
of  the  91 A  model  allows  11  more  stations  to  be 
considered.  If  the  comparisons  were  made  being  more 
selective  on  Doppler  stations  to  be  used  (at  least  3S 
passes  are  observed  after  June  1977)  the  mean 
differences  given  in  Table  1  become  more  positive  by 
IS  cm  and  the  standard  deviation  decreases  to  ±87  cm. 
Assuming,  somewhat  incorrectly,  that  the  undulation 
differences  arc  independent,  the  standard  deviation  of 
the  mean  difference  is  approximately  ±3  cm. 

Table  1 

Mean  and  Standard  Deviation  of  the  Geometric  and 
Gravimetric  Undulation  Difference  for  a  Global  Set  of 
Doppler  Stations 


GcQDQifintial 

Model 

OSII91A 

JGM.2/OSU91A 

Mean  Diff. 

•9  cm 

-15  cm 

SD  Diff. 

1100  cm 

99  cm 

No.  Stat. 
Accented 

1428 

1417 

Results  •  England/Germany 

The  next  step  was  to  examine  stations  that  could  be 
closely  associated  with  two  regions  with  different 
vertical  datums.  The  first  area  was  England  with 
stations  selected  in  the  region  50®  to  60®  N,  355®  to 
2®E.  We  assume  that  the  orthometric  heights  we  have 
are  given  in  the  Newiyn  Datum  (Ashkenazi  et  al., 
1990).  The  second  area  was  Germany  with  stations 
selected  in  the  region  48®  to  54®  N,  6®  to  13®E.  We 
assume  that  the  orthometric  heights  we  have  refer  to  the 

vertical  datum  described  in  Torge  (1980,  Section  6.2.3). 

1 

The  results  of  the  comparisons,  based  on  the  2m 
rejection  criteria,  arc  given  in  Table  2.  The  mean 
differences  imply  that  the  Newiyn  Datum  is  -87  cm 
below  the  geoid  while  the  German  datum  is  4  cm  above 
the  geoid.  The  difference  implies  that  the  German 
vertical  datum  reference  surface  is  91  cm  above  the 
Newiyn  Datum. 


Mean  and  Standard  Deviation  of  the  Geometric  and 
Gravimeuic  Undulation  Difference  for  Doppler 
Stations  in  England  and  Germany 


England 

Germapy 

Mean  Diff. 

-87  cm 

4  cm 

SD  Diff. 

±64  cm 

±75  cm 

No.  Stat. 

29 

28 

Accepted 

An  independent  estimate  of  the  above  can  be  obtained 
as  follows.  Willis  et  al  (1989)  estimate  (on  the  basis  of 
several  sources  of  information)  that  the  origin  of  the 
French  vertical  datum  (IGN69)  is  30  cm  ±8  cm  above 
the  origin  surface  of  the  Newiyn  Datum  (ODN).  (Dn  the 
basis  of  the  comparison  of  heights  near  the 
French/German  border  (Boucher  1993,  private 
communication),  the  German  system  origin  is  50  cm 
above  the  IGN69  system  (i.e.  H  (Germany)  =  H 
(IGN69)  -  50  cm).  This  numerical  value  can  be  caused 
not  only  by  vertical  datum  differences,  but  also 
differences  in  normal  heights  and  orthometric  heights. 
The  lauer  we  will  ignore  for  this  paper.  Taking  the  30 
and  50  cm  into  account  would  imply  that  the  German 
origin  surface  is  80  cm  above  the  Newiyn  Datum 
surface.  This  number  can  be  contrasted  with  the  91  cm 
obtained  through  the  Doppler  station  analysis. 

Results  in  the  United  States 

Ihe  orthometric  heights  for  the  Doppler  stations  in  the 
United  States  are  given  with  respect  to  the  National 
Geodetic  Vertical  Datum  of  1929  (NGDV29).  This 
datum  was  defined  through  a  general  adjustment  in 
1929  with  height  constraints  imposed  at  26  dde  gauge 
stations.  The  orthometric  heights  were  computed  using 
normal  corrections  since  no  gravity  observations  were 
available  along  the  leveling  lines. 

A  new  adjustment  of  leveling  data  has  led  to  the  Ntnth 
American  Vertical  Datum  of  1988  (Zilkoski,  Richards, 
and  Young,  1992).  This  datum  is  based  on  geopotential 
numbers  fixing  the  height  at  a  primary  tidal  benchmark. 
Although  we  did  not  have  heights  for  the  Doppler 
stations  directly  in  the  NAVD88  system,  a  datum 
conversion  software  program  (VERTCON,  version  1.0. 
Februap'  1993)  has  been  developed  at  the  National 
Geodetic  Survey.  This  software  was  used  to  convert  the 
NGVD29  heights  into  NAVD88  heights  using  a 
transformation  procedure  in  lieu  of  the  direct 
adjustment  of  the  points  in  NAVD88. 

The  stations  to  be  used  in  this  analysis  were  initially 
selected  in  the  region;  30®  to  45®N,  240®  to  290®E. 
This  station  set  was  then  thinned  to  eliminate  stations 
geographically  close  so  that  a  file  containing  350 


stations  was  analyzed.  The  results  of  the  undulation 
comparison  with  the  JGM-2/OSU91A  potential 
coefficient  model  and  a  2m  rejection  criteria  arc  given 
in  Table  3. 

Table  3 

Mean  and  Standard  Deviation  of  the  Geometric  and 
Gravimetric  Undulation  Difference  for  Doppler 
Stations  in  the  United  States 


Vertical  Datum 
NGVD29  NAVD88 
Mean  Diff.  -26  cm  -72  cm 

SDDiff.  ±83  cm  ±100  cm 

No.  Stations  321 _ 321 


The  mean  differences  for  NGVD29  and  NAVD88  differ 
by  46  cm.  Somewhat  surprising  is  the  somewhat  poorer 
fit  (100  cm  vs  83  cm)  when  the  NAVD88  datum  is 
used.  One  would  suspect  that  the  NAVD88  heights 
should  be  more  consistent  with  the  ellipsoidal 
hcighu/gravimctric  undulations  than  the  NGVD29 
heights  which  are  given  in  a  distorted  datum. 

Figure  4  of  (ibid.)  show  the  height  differences  between 
NAVD88  and  NGVD29,  This  map  shows  that  the 
heights  with  respect  to  the  two  datums  are  quite 
consistent  in  the  eastern  half  of  the  country.  Due  to 
this,  additional  tests  were  made  restricting  the  analysis 
to  stations  cast  of  2S4*’E.  The  results  of  this  analysis 
are  given  in  Table  4.  We  see  a  substantial  change  (from 
Table  3)  of  the  mean  difference  for  stations  on 
NAVD88  (-72  cm  (Table  3)  to  -38  cm  (Table  4)).  The 
corresponding  change  in  the  NGVD29  difference  is 
only  8  cm.  The  standard  deviation  of  the  undulation 
differences  in  going  from  NGVD29  to  NAVD88  is 
smaller  than  the  corresponding  case  in  Table  3. 

Table  4 

Mean  and  Standard  Deviation  of  the  Geometric  and 
Gravimetric  Undulation  Difference  for  Doppler 
Stations  in  the  Eastern  United  States 


YcnicaTPaLum 

NGVD2^ 

NAVD88 

Mean  Diff. 

-18cm 

-38  cm 

SD  Diff. 

±86  cm 

±92  cm 

No.  Stations 

180 

180 

These  computations  were  repeated  restricting  the 
analysis  to  stations  where  the  number  of  satellite  passes 
was  greater  than  35  and  the  observation  data  was  after 
June  15, 1978.  The  standard  deviations  would  typically 
be  reduced  by  about  10%.  The  mean  differences  would 
change  by  7  to  12  cm. 


The  results  of  Table  3  make  it  unclear  as  to  the 
appropriate  mean  difference  to  be  associated  with  the 
two  vertical  datums.  One  can  argue  that  a  single 
parameter  can  not  be  used  to  represent  the  difference 
between  the  NGVD29  reference  surface  and  the  geoid 
because  of  the  manner  in  which  NGVD29  was  adjusted. 
However,  there  is  only  a  8  cm  difference  between  the 
total  station  set  and  the  eastern  region  data.  On  the 
other  hand,  one  would  expect  the  greatest  consistency 
with  the  NAVD88  system.  Unfortunately  the  mean 
difference  differs  by  34  cm  from  the  total  station  set  in 
the  eastern  station  set.  With  no  other  information 
available,  one  perhaps  should  use  the  results  given  in 
Tabic  3  understanding  the  value  for  NAVD88  is  subject 
to  change  when  more  accurate  heights  in  the  NAVD88 
system  are  used  in  the  analysis.  Specifically  from  Table 
3  one  estimates  that  the  NGVD29  reference  surftce  is 
located  26  cm  below  the  geoid  while  the  NAVD88 
reference  surface  is  72  cm  below  the  geoid. 

Results  in  Australia 

The  heights  available  at  the  Doppler  stations  are 
approximate  orthometric  heights  referred  to  the 
Australian  Height  Datum  (AHD).  The  development  of 
this  datum  is  described  by  Leppert  (1974)  and  Rizos  et 
at.  (1991).  The  levelling  adjusunent  was  carried  out  in 
1971  using  height  constraints  at  30  tide  gauge  stations 
around  the  coastline  of  Australia. 

The  Doppler  stations  were  well  distributed  around 
Australia.  The  results  of  the  analysis  are  given  in  Table 
5  (AHD  (M)).  The  mean  difference  of  -68  cm  suggests 
that  the  average  AHD  reference  surface  is  located  68 
cm  below  the  geoid  as  defined  by  the  JGM-2/OSU91A 
potential  coefficient  model.  This  magnitude  conflicts 
with  estimates  of  dynamic  topography  which  suggests 
that  mean  sea  level  in  the  vicinity  of  Australia  is  dtove 
the  geoid.  (For  example,  sec  Figure  1 1  in  Rapp,  Wang, 
Pavlis(1991)). 

Table  5 

Mean  and  Standard  Deviation  of  the  Geometric  and 
Gravimetric  Undulation  Difference  for  Doppler 
Stations  in  Australia  and  Tasmania 


AHD  (M)  AHnm 
Mean  Diff.  -68  cm  -98  cm 
SD  Diff.  96  66 

No.  Stations _ 85 _ 4 


An  attempt  was  made  to  determine  the  difference 
between  the  AHD  as  used  in  the  mainland  and  the  AHD 
used  in  Tasmania.  This  test  was  a  follow-up  to  the 
study  described  by  Rizos,  Coleman  and  Ananga  (1991). 
The  first  step  was  to  use  four  Doppler  stations  in 
Tasmania  that  yielded  a  mean  difference  of  -98  cm  with 
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a  standard  deviation  of  ±66  cm  (see  Table  5  AHD(T)). 
This  result  implies  that  the  AHD  on  the  mainland  hiu 
us  reference  surface  30  cm  higher  than  the  AHD  in 
Tasmania.  This  value  is  similar  to  an  updated  (from 
Ri/os,  Coleman,  Ananga  (1991))  estimate  of  40  cm 
(Coleman,  1993,  pnvate  communication).  A  caveat 
lierc  is  that  this  analysis  uses  data  distributed  over  all  of 
Australia  while  the  Rizos,  Coleman,  Ananga  study 
considered  the  difference  between  the  AHD  (Victoria 
(.  oast)  and  AHD  (Tasmania). 

Conclusions 

This  paper  has  demonstrated  a  procedure  that  estimates 
the  separation  between  reference  surfaces  defined  by 
several  vertical  datums  and  the  geoid  that  is  defined  by 
geoid  undulations  computed  from  the  JGM-2  (degree  2 
to  70)  and  OSU91A  (degree  71  to  360)  potential 
coefficient  model.  The  ellipsoidal  heights  of  the 
stations  arc  determined  from  Doppler  positioning 
lochniques  after  a  uansformation  to  the  1TRF90  system. 
Although  the  accuracy  of  the  vertical  position  derived 
from  the  Doppler  positions  is  on  the  order 
of  ±80  cm,  the  large  number  of  stations  on  some 
s  crtical  datums  makes  these  tests  meaningful. 

In  the  analysis  wc  analyzed  data  that  can  be  associated 
with  vertical  datums  in  England,  Germany,  Australia, 
and  the  United  States.  We  assume  that  the  difference 
between  the  geoid  and  the  reference  surface  of  the 
vertical  datum  can  be  modeled  by  a  single  parameter. 
This  would  be  appropriate  if  the  datums  were  uniquely 
defined  by  a  reference  surface  through  a  specified  point 
However,  the  adjustment  of  several  of  the  datums  was 
carried  out  by  fixing  heights  at  tide  gauges  distributed 
around  the  borders  of  the  country.  This  fixing 
introduces  a  distortion  in  the  reference  surface 
(Laskowski,  1983)  with  respect  to  a  single  point  related 
The 

reference  surface.  Our  procedure  then  calculates  a 
separation  for  some  average  reference  surface 
representative  of  the  area  in  which  the  stations  are 
given. 

The  results  for  the  differences  in  the  origins  of  the 
vertical  datums  studied  for  this  paper  afiq  shown  in 
Figure  I.  Not  shown  are  the  results  for  AHD 
(Tasmania)  because  of  the  few  stations  available.  Two 
results  for  NAVD88  (U.S.)  are  shown  because  of  the 
large  differences  in  the  result.  One  estimate  is  based  on 
the  complete  station  set  after  conversion  from  NGVD29 
while  the  second  estimate  is  based  only  on  stations  in 
the  eastern  part  of  the  US.  From  Figure  1  one  clearly 
secs  the  differences  of  the  reference  surfaces  of  a  few 
vertical  datums.  It  should  be  emphasized  that  the 
numerical  value  shown  in  tables  and  Figure  1  are 
dependent  on  many  factors  (station  reference  frame, 
geoid  undulation  model,  equatorial  radius,  etc.). 


The  results  (Table  3/4)  for  the  two  vertical  datums  in 
the  U.S.  are  somewhat  troubling.  Why  should  the 
standard  deviation  of  fit  be  poorer  when  the  NAVD88 
(a  more  rigorous  datum)  heights  are  used  as  opposed  to 
the  NGVD29  values?  Why  should  the  mean  differences 
change  so  much  (44  cm)  when  the  results  for  the 
complete  station  set  vs  the  eastern  US  station  set  of 
NAVD88  are  considered?  Additional  study  with  more 
information  on  NAVD88  heights  is  needed. 

Another  area  of  concern  is  the  magnitude  of  the 
geoid/datum  separation  for  the  Australian  Height 
Datum.  On  the  basis  of  oceanographic  determinations 
of  mean  sea  level  one  would  expect  the  mean 
differences  to  be  somewhat  (-20  cm)  positive  instead  of 
the  -68  cm  found  in  these  results.  On  the  encouraging 
side  was  the  good  agreement  between  the 
AHD(mainland)/AHD(Tasmania)  difference  between 
the  30  cm  (AHD  (M)  higher))  and  the  recent  estimate 
by  Coleman  of  40  cm. 

To  resolve  some  of  these  questions,  more  accurate 
compulations  are  needed  that  would  use  more  precise 
geocenuie  positions  of  stations  on  various  vertical 
datums  and  more  refined  computation  for  geoid 
undulation  values  at  the  stations.  Procedures  to  use  this 
data  for  vertical  datum  definition  and  connections  are 
described  in  Rapp  and  Balasubramania  (1992). 

And  finally,  the  results  of  this  study  imply  that  height 
systems  can  be  inconsistent  by  one  or  two  meters.  This 
is  not  new  information  as  it  is  well  known  from  studies 
of  dynamic  height  variations  in  the  oceans.  There 
seems  to  be  a  clear  need  to  develop  a  unique  height 
system  that  can  be  tied  to  the  ideal  reference  surface, 
the  geoid.  Given  a  sufficiently  accurate  representation 
of  geoid  undulations,  orthometric  heights  can  be 
derived  from  an  ellipsoidal  height  computed  from  space 
positioning.  Additional  consideration  is  needed  to 
define  the  type  of  orthomctric  height  wanted,  or  to 
implement  the  determination  of  a  normal  height  through 
a  height  anomaly  calculation  instead  of  a  geoid 
undulation  computation.  Whatever  the  case,  we  are 
approaching,  although  not  yet  there,  the  situation  of 
having  suffleient  gravity  field  information  and 
ellipsoidal  height  values  to  ^Hne  a  world  height  system 
at  the  50  -  100  cm  accuracy  level. 
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Fig.  1.  Summary  of  Reference  Surface  Differences  with  Respect  to  the  Geoid.  Units  are  cm. 
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3.  Personnel 

The  principal  investigator  of  this  project  is  Richard  H.  Rapp,  Professor  Emeritus, 
Depanment  of  Geodetic  Science  and  Surveying.  Mr.  Nagarajan  Balasubramania,  Graduate 
Research  Associate,  also  panicipated  in  the  research  under  this  project.  Mr.  Balasubramania  is  a 
Ph.D.  candidate  at  The  Ohio  State  University. 

4.  Other  Activities 

During  this  project.  Dr.  M.  Kumar,  of  the  Defense  Mapping  Agency,  visited  us  to  learn 
about  the  research  activities  taking  place.  The  first  visit  was  on  October  22,  1992  while  the 
second  visit  was  on  March  23,  1993.  For  both  visits,  presentations  were  made  to  Dr.  Kumar. 
Copies  of  the  presentations  were  sent  to  Anestis  Romaides  at  Phillips  Laboratories  (Code  LW). 

5.  Travel 

In  August  1993,  Dr.  Rapp  traveled,  with  support  from  this  project,  to  the  NASA/Goddard 
Space  Flight  Center  for  a  meeting  with  NASA  and  DMA  representations.  This  meeting 
concerned  the  possible  joint  activities  of  NASA,  DMA,  and  Ohio  State,  in  the  production  of  a 
new  gravity  mt^el  that  could  be  used  for  precise  geoid  determination.  This  travel  was  the  only 
travel  supported  by  the  project.  A  follow-up  meeting  on  the  joint  gravity  field  venture  was  held 
in  October  1993  at  NASA/Goddard.  This  meeting  was  also  attended  by  Dr.  Rapp.  No  travel 
funds  from  this  project  were  used  for  this  travel. 

6.  Presentations 

Dr.  Kumar  and  Dr.  Rapp  submitted  an  abstract  for  a  presentation  at  the  1993  Fall  Meeting 
of  the  American  Geophysical  Union.  The  abstract,  as  it  appears  in  EOS,  Trans.  American 
Geophysical  Union,  Vol.  74,  No.  43,  October  26, 1993,  is  as  follows: 

Unification  of  World  Height  Systems 

M  Kumar  (Defense  Mapping  Agency,  Fairfax,  VA,  22031;  tel:  703-487-8151) 

R  H  Rapp  (Dept,  of  Geodetic  Science  and  Surveying,  The  Ohio  State  University,  Columbus,  OH 
43210;  tel:  614-292-6005;  e-mail;  RHRAPP@OHSTMVSA.ACS.OHIO-STATE.EDU) 

There  are  several  hundred  venical  datums  used  worldwide  and  in  some  areas  there  is  no 
information  available  that  can  be  used  to  determine  a  height  in  the  local  system.  Consequently, 
there  is  a  substantial  need  to  improve  our  vertical  positioning  capability  for  geodetic,  mapping, 
and  charting  purposes.  Numerous  theories  have  been  developed  that  enable  connection  and 
establishment  of  vertical  datums.  Such  theories  require  data  ^ven  in  a  specific  datum  with 
accurate  space  positioning  techniques  combined  with  dense  gravity  information.  An  alternative 
procedure  proposed  in  this  paper  is  to  introduce  the  geoid  as  the  unifying  global  vertical 
reference  surface. 

To  locate  the  geoid  with  respect  to  a  reference  ellipsoid,  geoid  undulations  must  be  accurately 
determined.  Today  a  global  geoid  undulation  accuracy  estimate  is  about  26  cm  in  ocean  areas. 
On  land,  the  achievable  accuracy  is  definitely  poorer  and  in  areas  with  scarce  gravity  data  the 
accuracy  may  be  even  2-5  m.  To  achieve  an  accurate  determination  of  the  global  geoid  for 
vertical  datum  purposes  improvement  in  our  knowledge  of  the  Earth's  gravity  field  is  needed. 

Thus,  as  a  first  step,  efforts  are  being  made  to  develop  a  new  degree  360  potential  coefficient 
model  combining  dl  the  available  satellite  and  terrestrial  gravity  data.  This  new  model  should 
then  enable  the  determination  of  geoid  undulations  to  an  accuracy  of  50  to  100  cm  in  land  areas 
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with  substantially  better  accuracy  over  ocean  areas.  The  new  geoid  is  then  intended  to  be 
implemented  as  the  vertical  reference  surface  for  the  World  Geodetic  System  (WGS)  1984. 
These  geoidal  undulations  would  also  enable  the  determination  of  orthometric  heights  from 
space  positioning  systems,  e.g.,  GPS.  with  respect  to  a  consistent  zero  reference  surface  globally. 
In  effect,  the  improved  geoid  will  result  in  the  realization  of  WGS84  as  a  complete  three- 
dimensional  reference  system.  The  connection  of  this  global  height  system  to  local  vertical 
datums  can  also  be  then  easily  achieved.  The  technical  approach  to  develop  such  a  united  world 
height  system  will  be  described  in  this  presentation. 

7.  Future  Activities 

Most  of  the  research  goals  defined  in  the  original  proposal  have  been  accomplished 
during  this  contract.  This  report  has  described  two  aspects  of  the  research  carried  out. 
Additional  details  of  the  research  by  Balasubramania  will  be  given  in  his  dissertation  which 
should  be  completed  by  the  2nd  quarter  of  1994. 

The  results  of  this  research  clearly  shows  the  need  to  more  clearly  define  our  vertical 
datums.  The  numerical  results  demonstrate  the  origin  inconsistencies  that  exist  in  a  number  of 
current  datums.  Although  these  inconsistencies  are  at  the  ±1  m  level,  they  are  of  increasing 
concern  in  today's  high  precision  geodetic  positioning.  We  have  seen  that  one  of  the  difficulties 
in  connecting  our  datums  in  the  lack  of  a  sufficiently  accurate  gravity  model  for  the  Earth. 
Current  models  can  be  used  to  derive  geoid  undulations  that  have  an  accuracy  of  ±0.5  to  ±5  m. 
With  improved  data  both  satellite  and  surface  gravity  data,  an  improvement  in  our  representation 
t)f  the  Earth's  gravitational  potential  can  be  anticipated.  This  improvement  could  lead  to  the  ideal 
reference  surface,  the  geoid,  to  be  used  for  a  global  vertical  reference  system.  Ultimately,  this 
surface  could  be  the  reference  for  a  World  Height  System. 
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